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We theoretically examine generation of spin-polarized current using multi-terminated quantum dot with spin- 
orbit interaction. First, a two-level quantum dot is analyzed as a minimal model, which is connected to N 
(> 2) external leads via tunnel barriers. When an unpolarized current is injected to the quantum dot from a 
lead, a polarized current is ejected to others, similarly to the spin Hall effect. In the absence of magnetic field, 
the generation of spin-polarized current requires N > 3. The polarization is markedly enhanced by resonant 
tunneling when the level spacing in the quantum dot is smaller than the level broadening due to the tunnel 
coupling to the leads. In a weak magnetic field, the orbital magnetization creates a spin-polarized current even 
in the two-terminal geometry (N = 2). The numerical study for generalized situations confirms our analytical 
result using the two-level model. 

PACS numbers: 72.25.Dc,71.70.Ej,73.23.-b,85.75.-d 



I. INTRODUCTION 

The spin-orbit (SO) interaction in semiconductors has been 
studied extensively from viewpoints of its fundamental re- 
search and application to spin-based electronics, "spintron- 
ics."- For conduction electrons in direct-gap semiconductors, 
an external potential U (r) results in the Rashba SO interac- 
tion^ 

H RS0 = i<r-[pxVU(r)], (1) 

where p is the momentum operator and er is the Pauli matri- 
ces indicating the electron spin s - er/2. The coupling con- 
stant A is markedly enhanced by the band effect, particularly 
in narrow-gap semiconductors, such as InAs and InSb.i£ The 
bulk inversion symmetry is broken in compound semiconduc- 
tors, which gives rise to another type of SO interaction, the 
Dresselhaus SO interaction. 6 It is given by 

H D so = -fi[px(Py ~ Pz)°"* + Py(P 2 z ~ P 2 x)°~y 

+Pz(P 2 x -Py)cr t ]. (2) 

In the presence of SO interaction, the spin Hall effect (SHE) 
is one of the most important phenomena for the application to 
the spintronics. It produces a spin current traverse to an elec- 
tric field applied by the bias voltage. There are two types of 
SHE. One is an intrinsic SHE, which is induced by the drift 
motion of carriers in the SO-split band structures. It creates a 
dissipationless spin current. 7 9 The other is an extrinsic SHE 
caused by the spin-dependent scattering of electrons by impu- 
rities. 10 Kato et al. observed the spin accumulation at sample 
edges traverse to the current, 11 which is ascribable to the ex- 
trinsic SHE with U (r) being the screened Coulomb potential 
by charged impurities in Eq. (HJ. 12 The extrinsic SHE is usu- 
ally understood semi-classically in terms of skew scattering 
and side-jump effect. 

In our previous studies, 1314 we theoretically examined the 
extrinsic SHE in semiconductor heterostructures due to the 
scattering by single artificial potential. The potential created 



by antidots, STM tips, and others, is electrically tunable. We 
adopted the quantum mechanical scattering theory for this 
problem. When the potential is axially symmetric in two di- 
mensions, U(r) with r = ^x 2 + y 2 in the xy plane, electrons 
feel the potential 

U eS = U{r) + U x {r)l z a z (3) 

in the presence of Rashba SO interaction. U\{r) = -AU'(r)/r 
has the same sign as U(r) if \U(r)\ is a monotonically decreas- 
ing function of r and A > 0. For electrons with cr z = 1, 
U e g = U(r) + U] (r)l z and as a result, the scattering for com- 
ponents of l z > (l z < 0) is enhanced (suppressed) by the SO 
interaction. For electrons with <x- = -1, the effect is oppo- 
site. This is the origin of the extrinsic SHE in two-dimensional 
electron system. We showed that the SHE is significantly en- 
hanced by the resonant scattering when U (r) is attractive and 
properly tuned. We proposed a three-terminal spin-filter in- 
cluding a single antidot. 

In the present study, we examine an enhancement of the 
"extrinsic SHE" by resonant tunneling through a quantum dot 
(QD) in multi-terminal geometries. The QD is a well-known 
device showing a Coulomb oscillation when the electrostatic 
potential is tuned by a gate voltageJ^ The number of elec- 
trons is almost fixed by the Coulomb blockade between the 
current peaks of the oscillation. At the current peaks, the res- 
onant tunneling takes place through discrete energy levels in 
the QD at low temperatures of k% T <k T with level broaden- 
ing T due to the tunnel coupling to the leads. Recently, the 
SO interaction in QDs of narrow-gap semiconductors and re- 
lated phenomena have been investigated intensively^^ 7 - We 
consider a situation in which a QD with SO interaction is con- 
nected to N (> 2) external leads via tunnel barriers. We use 
the term SHE in the following meaning: When an unpolarized 
current is injected to the QD from a lead (lead S), polarized 
currents are ejected to the other leads [Dl,. . .,D{N - 1)]. In 
other words, the QD works as a spin filter. We assume that 
the SO interaction is present only in the QD and that the av- 
erage of level spacing in the QD is comparable to the level 
broadening F (~ 1 meV), in accordance with experimental 
situations. 16 Thus the transport takes place through single or 
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a few energy levels in the QD around the Fermi level Bp in 
the leads. The strength of SO interaction Aso [absolute value 
of /iso in Eq. ©] is approximately 0.1 ~ 0.2 meV for InAs 
QDsZM and 0.23 meV for InSb QDs£- 

Our purpose is to elucidate the mechanism of SHE at a QD 
with discrete energy levels. Consider an electron with spin- 
up or -down injected to the QD from lead Dl (electric current 
flows from the QD to lead Dl). The SO interaction in the 
QD mixes a few energy levels around sp in a spin-dependent 
way [a rotation in the pseudo-spin space of the levels; see Eq. 
10], whereas the tunnel coupling to lead D2 mixes the levels 
differently in a spin-independent way. The interference be- 
tween the mixings results in the spin-polarized electron going 
out to lead S. To simply clarify the spin-dependent transport 
processes, we neglect the electron-electron interaction. We 
focus on the current peaks of the Coulomb oscillation where 
the interaction is not qualitatively important. 

First, we examine a two-level QD as a minimal model and 
present an analytical expression for the spin-dependent con- 
ductance. We assume single conduction channel in each of N 
leads. In the absence of magnetic field, we show that three or 
more leads (N > 3) are required to generate the spin-polarized 
current. We observe a large spin polarization by the resonant 
tunneling at the current peak when the spacing A between the 
two levels in the QD is smaller than T. Although the SHE 
at a QD seems quite different from the SHE by an impurity 
potential, the condition of A < F would correspond to the de- 
generacy for the virtual bound states with +L [see Eq. (0)]. 
The preliminary results of this part in the present paper were 
published in our previous paper.— 

Second, we analyze the transport through the two-level 
QD in a weak magnetic field. The orbital magnetization 
is taken into account to the first order of magnetic field, 
whereas the Zeeman effect is neglected. We find the creation 
of spin-polarized current in a conventional geometry of two- 
terminated QD (N = 2) with finite magnetization b [see Eq. 
(0; b ~ Hu> c with cyclotron frequency u> c = \e\B/m*] and en- 
hancement of the polarization when \b\ is comparable to the 
strength of the SO interaction Aso (magnetic field of B ~ 40 
mT). This is ascribable to the interference between the spin- 
dependent mixing of energy levels in the QD by the SO inter- 
action and spin-independent one by the orbital magnetization. 

Finally, our analytical results for the two-level QD are con- 
firmed by numerical study on the QD with several energy lev- 
els. A QD with tunnel barriers to N leads is modeled on a two- 
dimensional tight-binding model. We observe spin-polarized 
currents for N = 3 (N = 2) in the absence (presence) of mag- 
netic field. The spin polarization is markedly enhanced at the 
current peaks when a few energy levels are close to each other 
around sp. 

We make some comments here, (i) Previous theoretical 
papers 29 32 concerned the spin-current generation in a meso- 
scopic region, or an open QD with no tunnel barriers, in which 
many energy levels in the QD participate in the transport. 
Since we are interested in the resonant tunneling through one 
or two discrete levels in the QD, our situation is different from 
that in the papers. 

(ii) The present work indicates a QD spin filter in multi- 



terminal (two-terminal) geometries without (with) magnetic 
field although we emphasize the fundamental aspect of the 
mechanism for the SHE at a QD. Note that our spin filter 
works only at low temperatures since the SHE stems from the 
coherent transport processes through the QD. Other spin filters 
were proposed using semiconductor nanostructures with SO 
interaction, e.g., three- or four-terminal devices related to the 
SHE j 13 i 14 ' 33 ~ — a triple-barrier tunnel diode, 39 quantum point 
contact, 40,41 and a three-terminal device for the Stern-Gerlach 
experiment using a nonuniform SO interaction^ 

(iii) We do not consider the electron-electron interaction 
in the present paper focusing on the current peaks of the 
Coulomb oscillation. In the Coulomb blockade regimes 
between the current peaks, the electron-electron interaction 
plays a crucial role. We examined the many-body resonance 
induced by the Kondo effect in the blockade regime with spin 
1/2 in the multi-terminated QD. We showed the generation of 
largely polarized current in the presence of the SU(4) Kondo 
effect when the level spacing is less than the Kondo temper- 
ature. 28 We also mention that an enhancement of SHE by 
the resonant scattering or Kondo resonance was examined for 
metallic systems with magnetic impuritiesi^ 3 ^ 

The organization of the present paper is as follows. In Sec. 
II, we explain a model of two-level QD connected to N exter- 
nal leads. Section III presents the analytical expressions for 
the spin-dependent conductance using the model of two-level 
QD. In Sec. IV, we study a generalized situation in which a 
QD with many energy levels is connected to N leads through 
tunnel barriers. We make a two-dimensional tight-binding 
model to describe the situation and perform a numerical study. 
The last section (Sec. V) is devoted to the conclusions and dis- 
cussion. 



II. MODEL OF TWO-LEVEL QUANTUM DOT 

In this section, we explain our model depicted in Fig. 1(a), 
in which a two-level QD is connected to N external leads. 

We start from a QD with SO interaction and magnetic field 
B in general. The electronic state in the QD is described by 
the Hamiltonian 



H 



<(>> 



(P - eA) 1 
2m* 



+ U(r) + H so (B) 



A_ + U(r)+ ^. B . l + Hso , 
Zm Zm 



(4) 
(5) 



where U(r) is the confining potential of the QD, m* is the 
effective mass of conduction electrons (m*/mo = 0.024 in 
InAs with mo being the electron mass in the vacuum), and 
A = (B x r)/2 is the vector potential. Assuming a weak mag- 
netic field, we neglect the term of A 2 and Zeeman effect. For 
the SO interaction, //so can be the Rashba and/or Dresselhaus 
interactions in Eqs. ([TJ and ©. Although p in //so should 
be replaced by (p - eA) in the presence of magnetic field, the 
terms of A in Hso(B) can be disregarded in the case of weak 
magnetic field (see Appendix A). 

The eigenenergies of p 2 /(2m*) + U(r) form a set of discrete 
energy levels {e,}. We examine the situation in which two en- 
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FIG. 1: Models for a quantum dot (QD) connected to N (> 2) leads. When an unpolarized current is injected to the QD from lead S, polarized 
currents are ejected to leads Dl to D(N - 1). (a) A QD with two energy levels. The tunnel coupling between level j (= 1, 2) in the QD and 
lead a is denoted by V a j. (b) A QD (shaded square region of W x W in area) connected to three leads (quantum wires of W in width) via 
quantum point contacts, (c) A QD connected to two leads via quantum point contacts in the presence of magnetic field. Models in (b) and (c) 
are represented on a tight-binding model by discretizing the two-dimensional space (xy plane). 



ergy levels, ej and s<i, are relevant to the transport. The other 
levels are located so far from the two levels that the mixing by 
//so or \e\hB ■ I / '(2m*) can be neglected. 

The wavefunctions of the states, (r|l) and (r|2), can be real 
since they are eigenstates of real operator, p 2 /(2m*) + U(r). 
Since the orbital part in //so is a pure imaginary operator, it 
has off-diagonal elements only; 



<2|ffso|l> = i*so-o72 



(6) 



with h S o = h R so + h DS0 . ih RS0 = (2A/H)(2\(p x 
V£/)|l) in the case of Rashba interaction, whereas i/iDSO,.i = 
(2A' IK){2\p x (pl - p 2 )\l), etc., in the case of Dresselhaus in- 
teraction. For the same reason, ■ - (2\B ■ l\2) = 
and 



l ^(2\B-l\l)=ib/2. 
m 



(7) 



We estimate the value of \b\ to be \e\HB/m* = ttu c , where 
o) c - \e\B/m* is the cyclotron frequency. When hco c = 0.2 
meV (- Aso), the corresponding magnetic field is B = 40 mT 
in the case of InAs. If the quantization axis of spin is taken in 
the direction of hso, the Hamiltonian in the QD reads 



<T= + \ 



b + crAso 
T 7 + T, 



, (8) 



where a V. and dj a are the creation and annihilation operators 
of an electron with orbital /' and spin <x, respectively, e — (e\ + 
E2)/2, A = £2 - £i, and Aso = l^sol- The Pauli matrices, r y 
and r,, are introduced for the pseudo-spin representing level 1 
or 2 in the QD. Note that the Hamiltonian in Eq. (O yields the 
energy levels of s+ -y/A 2 + (b + cr Aso) 2 /2 for cr = + or - in an 
isolated QD; the Kramers degeneracy holds only with b = 0. 
Although the average of level spacing in a QD is assumed to 
be 6 ~ 1 me V, the spacing A between a specific pair of levels 
fluctuates around 6. A is fixed while the electrostatic potential, 
and hence the mean energy level e, is changed by tuning the 
gate voltage. 



The state | j) in the QD is connected to lead a by tunnel cou- 
pling, V a j (j = 1,2), which is real. The tunnel Hamiltonian 

is 



2 N 



/= 1 a= 1 k,o~ 
N 

= V a[{ e a,idl a + e aa d l 2ir )a aLo - + h.c.j, (9) 

where a a k,<7 annihilates an electron with state k and spin cr in 
lead a. V a = V(Vi) 2 + (V a , 2 ) 2 and e aJ = V aJ /V a . We in- 
troduce a unit vector, e a = (e j,c ff ,2) T . V a is controllable by 
electrically tuning the tunnel barrier, whereas e a is determined 
by the wavefunctions (r| 1 ) and (r|2) in the QD and hardly con- 
trollable for a given current peak. It should be mentioned that 
{e a } and A vary from peak to peak in the Coulomb oscilla- 
tion. We can choose a peak with appropriate parameters for 
the SHE in experiments. 

We assume a single channel of conduction electrons in the 
leads. The total Hamiltonian is 



EkC ak,cr Cak ' CT + ^ dot + 



(10) 



k,t 



The strengths of tunnel couplings to lead a are character- 
ized by the level broadening, F a = 7iv a (V a ) 2 , where v a is the 
density of states in the lead. We also introduce a matrix of 
f = E<r f with 



(ID 



An unpolarized current is injected into the QD from a 
source lead {a =S) and output to other leads [Dn; n = 
1, • • • ,(N — 1)]. The electrochemical potential for electrons 
in lead S is lower than that in the other leads by |e|Vbias- The 
transport through a QD in the multi-terminal geometry can be 
formulated following the paper by Meir and Wingreen^ just 
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FIG. 2: Spin-dependent conductance Gi,+ in the model of two-level quantum dot in the three-terminal geometry, as a function of mean energy 
level, e = (£!+ No magnetic field is applied. Solid (broken) lines indicate the conductance Gi + (G t ) for spin cr = +1 (-1) in the 

direction of hso (see Sec. II). The level spacing in the quantum dot is A = £2 - £1 = 0.21" (left panels) and T (right panels). The level broadening 
by the tunnel coupling to leads S and Dl is T s = T D1 = T (es,i/ e s,2 = 1> e Du/ e Di,2 = _ 2/3), whereas that to lead D2 is (a) r D2 = O.ir, (b) 
0.5r, (c) T, and (d) 2T (eD2,i/fiD2,2 = 2). The strength of spin-orbit interaction is fixed at Aso = 0.2T. 



as in the two-terminal geometry. The current with spin cr - + 
from lead a to the QD is written as 

W = ^ j deTr { f - [/«( £ ) (<% " G l) + > ( 12 ) 

where G^, GjJ., and Gf T are the retarded, advanced, and lesser 
Green functions in the QD, respectively, in 2x2 matrix form in 
the pseudo-spin space. f a (s) is the Fermi distribution function 
in lead a. 

Although the current formula in Eq. ( TTZl ) is applicable in 
the presence of electron-electron interaction in the QD, it is 
simplified in its absence. Then, & a - G'% = -liG^TG^ and 
G^j = 2iGJ r (2 ff r»/ff)G|J-. The substitution of these relations 
into Eq. ( TTZb yields 

W = jj ds L/b(«) - /s( £ )] Tr (<%f D„c; r f s ) , 

where /d«( s ) = /d(s)- At T = 0, the conductance into lead 
Dn with spin cr is given by 



G„ jC r 



d/n 



dV bi; 



= ^Tr[G^(e F )f D „G^( eF )fs], (13) 



V biils =0 



where the QD Green function is 



s - s d + f 

■ b±A s , 



■ b±A s , 



+ ir 



(14) 



HI. ANALYTICAL RESULTS 

We analyze the model of two-level QD, introduced in the 
previous section. We show analytical expressions for the spin- 
dependent conductance in the absence and presence of mag- 
netic field, respectively. 



A. In absence of magnetic field 

We begin with the case of b — 0, or in the absence of mag- 
netic field. From Eqs. ( fT3l and (1141 . we obtain 



Go cr — 



e 2 4F s r D „ r^ ( 
h \D\ 2 



£f - £ - — |eD«,tes,t 



,(2) 



-i js F - e + — |e D ,7.2''s.: 
A S0 , , 



+ 2^ F„(e D „ x e„) z (e s x e ff )- 



(15) 



(16) 



2 



(17) 
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FIG. 3: Spin-polarized conductance, Gi_ + - G[__, in the model of 
two-level quantum dot in the three-terminal geometry, as a function 
of e. The level spacing in the quantum dot is (a) A = £2 - £1 = 0.2r 
and (b) T. T D2 = 0.1X (solid line), T (broken line), and 2T (dotted 
line). The other parameters are the same as in Fig. 2. 

where D is the determinant of [G^sf)] 1 in Eq. (14\ , which 
is independent of cr. (a x b) z = a\b2 - a2b\. 

Let us consider two simple cases. (I) When A » T a and 
Aso, G«,<r consists of two Lorentzian peaks as a function of 
s, reflecting the resonant tunneling through one of the energy 
levels, £i_2 = e + A/2: 

G„,_-r s r D „^ ( ^_ gF)2 + (r ^ )2 . (is) 

Here, = 2„7rV(i(Kj) 2 is the broadening of level j (jj 
component of matrix F). In this case, the spin-polarized cur- 
rent [oc (G„ j+ - G„-)] is very small. A should be comparable 
to or smaller than the level broadening to observe a consid- 
erable spin current. (II) In a two-terminated QD (iV = 2), 
the second term in g^ ± vanishes. Since g^\ = g^_, no spin- 
polarized current is generated: 49 Three or more leads are re- 
quired to generate a spin-polarized current, as pointed out by 
other groups i 29 ' 31 ' 50 

We examine G\ t± in the three-terminated system (N = 3) in 

(2) 

the rest of this subsection. Thengj^ = [+(Aso/2)(es xe D\)z + 
Td2(^di x^Di)z( e s x^D2)z] 2 - We exclude specific situations in 
which two out of es, e^u and eoi are parallel to each other. 
The conditions for a largely spin-polarized current are as fol- 
lows: (i) A < (level broadening), as mentioned above. Two 
levels in the QD should participate in the transport, (ii) The 
Fermi level in the leads is close to the energy levels in the QD, 
ep ~ s (resonant condition), (iii) The level broadening by the 
tunnel coupling to lead D2, Fd2, is comparable to the strength 
of SO interaction Aso- 

Figures 2 and 4 show two typical results of the conductance 
Gi. ± as a function of s. In g^\ eDi,i<?s.i and eDi,2<?s,2 have 
different (same) signs in Fig. 2 (Fig. 4). Therefore, gf> = 
has no solution (a solution) in -A/2 < s - ep < A/2. 

In Fig. 2, the conductance shows a single peak. We set 
Fs = Tdi = T and change F^i from (a) 0.1T to (d) 2r. When 
A = 0.2F (left panels), we observe a large spin polarization 
around the current peak, which clearly indicates an enhance- 
ment of the SHE by the resonant tunneling [conditions (i) and 
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FIG. 4: Spin-dependent conductance Gi -± in the model of two-level 
quantum dot in the three-terminal geometry, as a function of mean 
energy level, s = (ei + No magnetic field is applied. Solid 

(broken) lines indicate the conductance Gi_ + (Gi _) for spin a = +1 
(-1) in the direction of hso (see Sec. II). The level spacing in the 
quantum dot is A = e 2 - £[ = 0.5T. The level broadening by the 
tunnel coupling to leads S and Dl is Ts = Tdi = T (es,i/«s,2 = 1> 
e D11 /e D i 2 = 2/3), whereas that to lead D2 is (a) r D2 = O.ir, (b) 
0.5r, (c) T, and (d) 2r (eD2,i/£D2,2 = 2). The strength of spin-orbit 
interaction is fixed at A so = 0.2r. Inset: Spin-polarized conductance, 
oc G+ - G_, as a function of s. T02 = OAT (solid line), 0.5r (broken 
line), and 2r (dotted line). 



(ii) ]. With increasing Fd2, the spin current increases first, 
takes a maximum in panel (c), and then decreases [condition 

(iii) ]. This means that the SHE is tunable by changing the 
tunnel coupling. When A = F (right panels), the SHE is less 
effective; spin polarization of P — (G\ t + - G\ _)/(Gi,+ + G\ _) 
around the current peak is smaller than in the case of A = 0.2r. 
However, a value of spin-polarized conductance, G n j+ - G„ _, 
is still large, as depicted in Fig. 3. 

In Fig. 4, the conductance G\ ± shows a dip at s * £f for 
small Fd2- The conductance dip is caused by the destructive 
interference between propagating waves through two orbitals 
in the QD. In the two-terminated QD without SO interaction 
(Fp2 = Aso = 0), the conductance Gi itr oc gy would com- 
pletely vanish at the dip, where the "phase lapse" of the trans- 
mission phase takes place. 51 As seen in Fig. 4, the conduc- 
tance dip changes to a peak with increasing Fd2 in the three- 
terminated QD. The SO interaction makes a large difference 
between Gi + and G\ - around the dip or peak, similarly to the 
case in Fig. 2. The spin-polarized conductance, Gi + - Gi _, 
shows a large peak there, as seen in the inset in Fig. 4. In 
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FIG. 5: Spin-dependent conductance Gi_ ± in the model of two-level quantum dot in the two-terminal geometry, as a function of mean energy 
level, £ = (£[+ £2)/2, in the presence of magnetic field. Solid (broken) lines indicate the conductance G\ t+ (Gi _) for spin cr = +1 (-1) in 
the direction of h so (see Sec. II). The level spacing in the quantum dot is A = e 2 - S\ = 0.2T (left panels) and T (right panels). The level 
broadening by the tunnel coupling to leads S and Dl is T s = r D [ = T (e & i /e S2 = 1, e ol ,t/e ol2 = -2/3). The orbital magnetization is (a) 
b = 0.021", (b) 0. IT, (c) 0.5r, and (d) T. The strength of spin-orbit interaction is fixed at A so = 0.2T. 



Fig. 4(a) with Fd2 = 0.1F, we find that the spin polarization 
of P = (Gi,+ - G\ _)/(Gi,+ + G\ _) is close to unity around the 
dip since G\- is almost zero. 



B. In presence of magnetic field 

Now we discuss the case with magnetic field: b + 0. The 
conductance into lead D« with spin cr is 



Gn o- — 



e 2 4F s F D n [ ;|i 
h \D a \ 2 



(19) 



where is the same as that in Eq. ([ToT i. whereas 



(2) 
?n,± 



b± A 



SO 



(e s x e D „) z 



+ ^ r a (e D „ x e a ) z (e s x e a ), 



(20) 



The determinant of [GJ,-(£f)] , D a , depends on cr in this case. 

In contrast to the case of b = 0, we observe the spin- 
dependent transport in a conventional geometry of two- 
terminated QD (N = 2). Then gf' ± = (b + A so ) 2 (e s x e m ) 2 /4. 
We expect a large spin polarization when (iii') b and Aso are 



comparable to each other, besides conditions (i) and (ii) in the 
previous sebsection are satisfied. 

We focus on the two-terminated QD (Af = 2) in this subsec- 
tion. Figures 5 and 7 exhibit the spin-dependent conductance 
Gj i± as a function of s. eDi,tes,t and eo\,2es,2 have different 
(same) signs in Fig. 5 (Fig. 7). We set Fs = Fqi = F, whereas 
the orbital magnetization is gradually increased from (a) to 
(d), or from (a) to (c). 

In Fig. 5, the level spacing in the QD is A = 0.2r in the 
left panels and A = T in the right panels. In the absence of 
magnetic field (b = 0), we do not observe the spin-polarized 
current in the two-terminal geometry, as discussed in the pre- 
vious subsection. With an increase in b, the difference be- 
tween Gi i+ and G\- increases, becomes maximal at b ~ Aso, 
and decreases [condition (iii')]. The SHE is more prominent 
for A = 0.2F than for A = F; the polarization P is larger in the 
former. 

Figure 6 shows the spin-polarized conductance, G\ )+ -G\-, 
as a function of s. We observe a large value even in the case 
of A = r if both the magnetic field and s are properly tuned. 

In Fig. 7, we observe a dip of conductance at e ~ sp. 
Around the dip, the spin-polarized current is largely enhanced 
as shown in the inset. In Fig. 7(a) with b = 0.1F, the spin po- 
larization of P — (Gi >+ - Gi _)/(Gi ;+ + Gi _) is close to unity 
because G\ _ almost vanishes. 
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FIG. 6: Spin-polarized conductance, Gi_ + - Gi__, in the model of 
two-level quantum dot in the two-terminal geometry, as a function of 
£, in the presence of magnetic field. The level spacing in the quantum 
dot is (a) A = e 2 - S\ = 0.2T and (b) T. The orbital magnetization b 
is b = 0.02r (solid line), O.ir (broken line), and 0.5r (dotted line). 
The other parameters are the same as in Fig. 5. 
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IV. NUMERICAL STUDY 

In the previous section, we have presented the analytical 
expressions for the spin-dependent conductance for the model 
of two-level QD. We have illustrated the generation of spin- 
polarized current in three- and two-terminated geometries in 
the absence and presence of magnetic field, respectively. In 
this section, we perform numerical studies for the QD with 
many energy levels to confirm our analytical results. A QD 
with tunnel barriers to N leads (N = 2, 3) in Figs. 1(b) and (c) 
is modeled on the tight-binding model in the xy plane. 



FIG. 7: Spin-dependent conductance Gi -± in the model of two-level 
quantum dot in the two-terminal geometry, as a function of mean 
energy level, e = (e\ + £i)/2, in the presence of magnetic field. Solid 
(broken) lines indicate the conductance Gi_ + (Gi _) for spin <x = +1 
(-1) in the direction of hso (see Sec. II). The level spacing in the 
quantum dot is A = e 2 - Si = 0.5r. The level broadening by the 
tunnel coupling to leads S and Dl is Ts = Tdi = T (es,i/«s.2 = 1> 
£Di,i/ e Di,2 = 2/3). The orbital magnetization is (a) b = OAT, (b) 
0.5r, and (c) T. The strength of spin-orbit interaction is fixed at 
Aso = 0.2T. Inset: Spin-polarized conductance, oc G+ - G_, as a 
function of s. b = OAT (solid line), 0.5r (broken line), and T (dotted 
line). 



A. Model 

In Figs. 1(b) and (c), N leads connect to a QD via tunnel 
barriers. The N leads are represented by quantum wires of 
width W with hard-wall potential at the edges. The electro- 
static potential in the QD (shaded square region of W x W) is 
changed by eV g . 

The tunnel barriers are described by quantum point contacts 
(QPCs). Along a quantum wire in the x direction, the QPC is 
described by the potential 47 



U(x,y;Uo) 



Uo 
2 



1 + cos 



+ £f 



Z 



£qpc 

2 



Wqpc 



at -Lqpc < x < Lqpc, where 



y±(x) 



w 
T 



1 - cos 



^QPC 



(21) 



(22) 



and 9{t) is a step function [6 = 1 for t > 0, 6 = for t < 0]. 
Uo is the potential height of the saddle point of QPC, whereas 
Lqpc and Wqpc characterize the thickness and width of the 
QPC, respectively. In the QD (shaded square region), the QPC 
potential is modified to U(x,y; Uq - eV g ) + eV g . In Fig. 1(b), 



we cut off the QPC potential at the diagonal lines of the square 
to avoid the overlap of two QPC potentials. 

As for the SO interaction, we consider the Rashba interac- 
tion caused by the QPC potential in Eq. d2TI ). that is, 



A 

Hso = -rcr z 
n 



dU 



Py 



dU 
dx 



(23) 



We choose the z direction for the spin axis (hso II Z direction). 

In the two-terminal geometry of Fig. 1(c), we consider a 
magnetic field perpendicular to the xy plane only in the region 
surrounded by dotted line. We adopt the vector potential of 
A = (-By, 0, 0) for the orbital magnetization and neglect the 
Zeeman effect. 

We discretize the two-dimensional space with QPC poten- 
tials and obtain the tight-binding model. We numerically eval- 
uate the spin-dependent conductance, using the calculation 
method in Appendix B. 

We consider the following situation. The width of quan- 
tum wires is W — 100 nm. The lattice constant of the tight- 
binding model is a — W/30 (number of sites is M = 29 in 
width of the wires). For the SO interaction, the dimension- 
less coupling constant is A — A/(2a 2 ) = 0.05, which corre- 
sponds to A = 1.171 nm 2 in InAs. 4 The Fermi wavelength 
and Fermi energy in the leads are fixed at Ap - W/3 and 
£p/t — 2 - 2 cos(2na / Ap) =* 0.382, respectively. (There are 
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FIG. 8: Numerical results of spin-dependent conductance G\ y± in 
model of Fig. 1(b), as a function of electrostatic potential in the quan- 
tum dot, eV g . No magnetic field is applied. Solid (broken) lines in- 
dicate the conductance Gi_ + (Gi_) for spin a = +1 (-1) in the z 
direction. The height of QPC potential is f7o = 0.8ef at the connec- 
tion to leads S and Dl, whereas (a) U /s F = 1.1, (b) 0.9, (c) 0.8, and 
(d) 0.6 at the connection to lead D2. 



six conduction channels in each lead. However, single chan- 
nel is effectively coupled to the QD owing to the QPC po- 
tential between the QD and lead.) For the QPC potential, 
Lqpc = Wqpc = /If- Uu = 0.8sf at the connection to leads 
S and Dl, whereas Uq at the connection to lead D2 is changed 
from U{)/sp = 1.1 to 0.6 to tune the tunnel coupling Td2 in 
the three-terminal geometry of Fig. 1(b). In Fig. 1(c), the 
magnetic field is applied up to Ha> c /sp = 30 x 10~ 4 , which 
corresponds to B 34 mT. 



B. NUMERICAL RESULTS 

Figure 8 presents the spin-dependent conductance G\ tCr in 
Fig. 1(b) of three-terminated QD, in the absence of magnetic 
field, cr — + 1 indicates the z-component of electron spin. The 
conductance shows a peak structure as a function of electro- 
static potential in the QD, eV s , reflecting the resonant tunnel- 
ing through discrete energy levels in the QD. Although this 
is similar to the Coulomb oscillation, the peak-peak distance 
is underestimated because we neglect the electron-electron in- 
teraction. 
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FIG. 9 : Spin-polarized conductance, G \ _+ - G \ _ , as a function of e V g 
in model of Fig. 1(b). The height of QPC potential at the connection 
to lead D2 is U /e f = 1.1 (solid line), 0.9 (broken line), 0.8 (dotted 
line), and 0.6 (dotted broken line). The other parameters are the same 
as in Fig. 8. 



The average of the level spacing is larger than the level 
broadening in Fig. 8. Therefore, the difference between G\ )+ 
and G\ _ is usually small. We observe a large spin-polarized 
conductance, Gi j+ - Gi,_, around some conductance peaks 
where a few levels should be close to each other around sp- 
Look at the conductance around eV g /sp = 0.25. With increas- 
ing the tunnel coupling to lead D2 by decreasing the height of 
QPC potential Uq, the spin-polarized conductance increases, 
becomes maximal, and decreases. This is in accordance with 
the analytical result in section III. A although it is hard to eval- 
uate the level spacing and signs of tunnel coupling around e F . 
Figure 9 plots Gi j+ - G\ - as a function of eV g , which seems 
complicated probably due to the interference among three lev- 
els in the QD. 

Figure 10 shows the spin-dependent conductance Gi jtr in 
Fig. 1(c) of two-terminated QD, in the presence of magnetic 
field. A large spin-polarized conductance is obtained around 
eVg/eF = -0.21. The difference of G\ + - G\ - is changed 
with increasing magnetic field perpendicular to the QD. As 
shown in Fig. 11, the absolute value of spin-polarized con- 
ductance increases, becomes maximal, and decreases, in ac- 
cordance with the analytical result in section III.B. 



V. CONCLUSIONS AND DISCUSSION 

We have studied the mechanism of SHE at a QD with dis- 
crete energy levels in multi-terminal geometries. We have 
considered a QD with SO interaction connected to exter- 
nal leads via tunnel barriers. When an unpolarized current is 
injected to the QD from a lead, a polarized current is ejected 
to others. > 3 (N > 2) is required in the absence (presence) 
of magnetic field for the generation of spin-polarized current. 

First, we have obtained the analytical expressions for the 
spin-dependent conductance using a minimal model of two- 
level QD. The SHE is markedly enhanced by the resonant 
tunneling when the level spacing in the QD is smaller than 
the level broadening due to the tunnel coupling to the leads. 
In the absence of magnetic field, the spin polarization can be 
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FIG. 10: Numerical results of spin-dependent conductance Gi.+ in 
model of Fig. 1(c), as a function of electrostatic potential in the quan- 
tum dot, eV g , in the presence of magnetic field. Solid (broken) lines 
indicate the conductance Gi,+ (Gi _) for spin <x = +1 (-1) in the z 
direction. The magnetic field is (a) H(i) c /ej, = 2 x lfr 4 , (b) 10 x 10~ 4 , 
and (c) 30 x 10~ 4 . 



tuned by changing the tunnel coupling to the lead other than 
source and drain leads in a three-terminal geometry. A weak 
magnetic field can tune the spin polarization in a two-terminal 
geometry. 

Second, we have performed numerical studies on the tight- 
binding model representing a QD with tunnel barriers to N 
leads. We have observed a large spin-polarized conductance 
at some current peaks when a few energy levels in the QD are 
close to each other around ep- The numerical results are in 
accordance with our analysis of the minimal model of two- 
level QD. 

In our calculation, we have neglected the electron-electron 
interaction. Therefore, our theory is applicable only around 
the current peaks of the Coulomb oscillation, where the in- 
teraction is not qualitatively important. We have also ne- 
glected the Zeeman effect. In spite of a large g-factor in InAs 
(\g\ ~ 10),— r 22 i 25 i 26 the Zeeman effect is smaller than the or- 
bital magnetization by one order of magnitude for B ~ 40 
mT, as estimated in Appendix A. In the absence of SO inter- 
action, the Zeeman effect splits a spin-degenerate level in the 
QD, which could result in the spin-polarized current by the 
resonant tunneling through one of the spin-split levels. In our 
situation, however, the spin splitting is much smaller than the 
level broadening, and hence the spin polarization by the Zee- 
man effect is negligibly small. 

We discuss a possible observation of the SHE at a QD. 
Since the measurement of a spin-polarized current is usu- 
ally difficult, an alternative method is desirable. Hamaya et 
al. fabricated InAs QDs connected to ferromagnets,— If a 
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FIG. 11: Spin-polarized conductance, Gi,+ - Gi,_, as a function 
of eVg in model of Fig. 1(c) in the presence of magnetic field. The 
magnetic field is Sw c /e F = 2 x 10~ 4 (solid line), 10 x 10~ 4 (broken 
line), and 30 x 10~ 4 (dotted line). The other parameters are the same 
as in Fig. 10. 



ferromagnet is used only for lead S and a normal metal or 
semiconductor is used for the other leads, an "inverse SHE" 
takes place. The electric current to lead Dl is proportional to 
(1 + p cos 0)Gi i+ + (1 - p cos 6)G\-, where p is the polariza- 
tion in the ferromagnet and 9 is the angle between the mag- 
netization and h$o. Thus G\ + and G\ _ can be evaluated by 
measuring the electric current with rotating the magnetization 
of the ferromagnet. 

The QDs are highly tunable systems. We believe that the 
detailed study of the generation of spin-polarized current at 
the QDs would contribute to the deeper understanding of the 
SHE. 
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Appendix A: Approximation in weak magnetic field 

In the presence of magnetic field and SO interaction, the 
Hamiltonian in Eq. © is approximated to that in Eq. (O for 
the following reasons. We consider the situation in which the 
average of level spacing in the QD is 5 ~ h 2 /(m*d 2 ) =2 1 meV, 
the strength of SO interaction is Aso - 0.2 meV, and magnetic 
field of hoj c ~ Aso or smaller. Here, d is the one-dimensional 
size of the QD, effective mass m*/mo - 0.024 in InAs, and 
(jd c = \e\B/m* is the cyclotron frequency. 

We choose the gauge of A = (6 x r)/2 in the model of 
two-level QD. The first-order in A in Hamiltonian © gives 
rise to the matrix element in Eq. (7). \b\ ~ Hu> c , as denoted in 
Sec. II. The second order in A yields e 2 / (8m*){i\(B x r) 2 \ j) ~ 



10 



(eBd) 2 /m* = (haj c ) 2 /6, which is smaller than the first-order 
term by ftcj c 1 5 <k 1 . 

We also neglect the Zeeman effect, Hz - g^B-a-/!, where 
fiB = \e\h/(2mo) is the Bohr magneton. We estimate the ef- 
fect to be \g\/i B B/2 = (\g\/4)(m"/m )Tiat c . Since \g\ ~ 10 in 
InAsr^ " 22 ' 25,26 The Zeeman term is smaller than the orbital 
magnetization \b\ by one order of magnitude. 

In the presence of magnetic field, p in H$o is replaced by 
(p - eA). In the case of Rashba interaction, 



Hrso(B) = j[(p-eA)xVU] 

n 

= j(p x Vt/) - ^[(Bxr)x Vf/] . 
n 2n 



(Al) 



The matrix element of the first term in Eq. ( IA1I ) is estimated to 
be (A/H)\(2\p x Vf/|l)| ~ (A/d 2 )6, whereas that of the second 
term is to be |e|/l/(2/z)Ki|(B x r) x VU\j)\ ~ (|e|AB/^)<5. The 
latter is smaller than the former by ha> c /6 <K 1, and thus it is 
safely disregarded. In the case of Dresselhaus interaction, 

i/ 

H DSO (B) = — [(n y n x n y -n z n x n z )cr x 
n 

+(jr z n y 7T z - n x 7T y n x )o- y 

+ (n x 7I z 7T x — 7T y 7I z 7Ty)o~ z ], (A2) 

where n - p-eA. The matrix element of the terms without A 
[Eq. (f2|i] is estimated to be (A' ft 2 /cP) and that of the first order 
in A is to be (A'H/d)\e\B. Again, the latter is smaller than the 
former by ftu) c /d <K 1. 



Appendix B: Numerical calculation of tight-binding model 



with B - \e\Ba 2 /h. B is related to the cyclotron frequency by 
B = hu> c /(4nt). The transfer term in the y direction is given 
by 



Tj,i-jj+i-± - 1 + iA(U j+i/2,/+i/2 - U j- 1/2,/+ 1/2)- 



(B3) 



To randomize the discrete energy levels in the QD, we in- 
troduce a uniformly distributed on-site energy Wij in the range 
of -W m J2 < w u < W mn /2. We choose W rdn '= 0.2e F . We 
disregard the SO interaction induced by the random potential. 

The spin-dependent conductance is numerically evaluated 
in the following way. First, we define the channels in the leads 
outside of the dotted line, which are represented by the quan- 
tum wires of width W = (M + \)a. Consider a quantum wire 
in the x direction. There are M channels, M con d of which are 
conduction modes and (M-M con d) are evanescent modes. The 
wavefunction of conduction mode fi {/u = 1, 2, • • • , M con d) is 
written as 



<M.A0 = exp(i^a;')^(0 



2 I nul 

sin 



M+ 1 



M+ 1 



(B4) 
(B5) 



with I = 0, 1, 2, • • • , M. The wavenumber k M satisfies e^C^) = 
£f, where the dispersion relation is given by 

E Jk) = 4t- 2t cos ( — | - 2tcos( lea). (B6) 
\M + 1 / 

The band edge, s^k = 0), is located below ep for the con- 
duction modes. The wavefunction of evanescent mode // 
in = Mcond + !,••■ , M) is written as 



</f/iU> = exp(-^a;>^(0, 



(B7) 



In the model of Figs. 1(b) and (c), we discretize the xy 
plane with QPC potentials and obtain the two-dimensional 
tight-binding model of square lattice. 46 The lattice constant 
is denoted by a. For the region surrounded by dotted line, the 
Hamiltonian is given by 



H = 



f Z( 4+ ^') 

j,l,<r 
j,l,cr 



C j,l;o- C j' l ' lr 



j+l,l;o-Cjj. a .Cj+ij. (r + Tjj-jJ+i ;(r Cjj. .Cj,M;a- + 



h(B)) 



where a j is the distance from the QD along the lead. The band 
edge is located above ep and k m is determined from ^(i/c^) = 

Next, we introduce the retarded Green function Go-(e) for 
the inside region of dotted line in Figs. 1(b) and (c). Here, 
cr = +1 represents the z component of spin, which is a good 
quantum number in Hamiltonian (IB lb . It is defined by 



G«r(e) = 



el-'Hc 



»=S,D« 



(B8) 



where cj / cr and Cjj ;cr are creation and annihilation operators of 
an electron at site (J, I) with z-component of spin cr = + 1 , re- 
spectively. The transfer integral is t = h 2 /(2m*a 2 ). Ujj repre- 
sents the QPC potential and electrostatic potential in the QD, 
at site (j, I) in units of t. The transfer term in the x direction is 
given by 



TjJ-j+\,l-± - {l +iA(Uj+\/2J+\/2 ~ C r ;'+l/2,/-l/2)} 



\2nBI 



where A = A/ '(2a 2 ) is a dimensionless strength of SO inter- 
action and f7 i+ i/2,/+i/2 is the potential at the middle point be- 
tween the sites (J, I) and (j + 1,1 + 1). The magnetic field in 
the QD is taken into account by the Peierls phase factor, e l2nBl 



where < H cr is the matrix of Hamiltonian with spin cr 
the self-energy due to the tunnel coupling to lead a{- 
and given by 



Z a = -tTlUKU- l T a . 

U = (Ki,«2, ••■ ,Um) is an unitary matrix, with 
u,j = (w„(l), 1^(2), • • • ,m^(M)) t in Eq. dB5t. A = 
(B2) diag(/li, /I2, • • ■ ,Am), where A^ = exp(ifc A ,a) for conduction 
modes and A^ = exp(-/^fl) for evanescent modes. T a is a 
coupling matrix between the edge of lead a to the consider- 
ing region; T a (l, a/) = 1 if site a/ is connected to the site / 
(= 1,2, • ■ • , M) at the end of the lead, T a (l, ai) = other- 

,i,io 46 



=. E ff is 
S,Dn) 



(B9) 
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The conductance from lead S to Dl can be evaluated sep- at T — 0, where T a = i[Y. a — Y. a ]/2. 
arately for cr = ±1 of the z component of spin. The spin- 
dependent conductance is calculated using the formula 

G u± = ^Tr[r D1 G ± (s F )r s Gl( £F )] (BIO) 
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